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The probability density function (PDF) for critical wavefunction amplitudes is studied in the three- 
dimensional Anderson model. We present a formal expression between the PDF and the multifractal spectrum 
/(a) in which the role of finite-size corrections is properly analyzed. We show the non-gaussian nature and 
the existence of a symmetry relation in the PDF. From the PDF, we extract information about f{a) at criticality 
such as the presence of negative fractal dimensions and we comment on the possible existence of termination 
points. A PDF-based multifractal analysis is hence shown to be a valid alternative to the standard approach 
based on the scaling of general inverse participation ratios. 
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The fluctuations and correlations of wave amplitudes are 
of primary importance for the understanding of many classi- 
cal and quantum systems. This is arguably most pronounced 
in the physics of Anderson localization [Jl. Here, recent ad- 
vancements in theory jHSH], experiments in classical jstl 
and quantum waves la III l8|] as well as numerical meth- 
ods 10, [Toll have led to unprecedented insights into the na- 
ture of the localization-delocalization transition. In contrast 
to the weak- or strong-disorder limits where the description 
of nearly-extended or strongly localized states is well-known, 
e.g. from random matrix theory ifTHl . the intensity distribu- 
tion at the metal-insulator transition is more involved due to 
the multifractal nature of the states iH [13, 12]. The pos- 
sibility of carrying out a multifractal analysis directly from 
the raw statistics of intensities j'i/'p, i.e. the probability den- 
sity function (PDF), is especially interesting since their jiistri- 
butions can be measured experimentally in classical 15] and 
quantum 10. E, 15] experiments. The PDF at criticality is 
closely related to the multifractal spectrum /(a). However, 
the numerical relation between the PDF and f{a) has not 
been completely elucidated. In this Letter we show how to 
obtain the multifractal spectrum based on the PDF. The PDF- 
to-/(Q;) connection is a numerically much sim pler procedure 
than the usual scaling of g-moments of 11611 . Further- 
more, it yields direct understanding of physical properties at 
criticality, such as the existence of a symmetry relation, the 
observation of negative fractal dimensions and the physical 
meaning of the possible appearance of termination points. We 
apply the PDF-based approach to the three-dimensional An- 
derson model within the Gaussian orthogonal ensemble, us- 
ing a large number of critical states at E = and very large 
system-sizes up to = 240^ l!^. 

At criticality, the intensity distribution has the scaling form 

VlM^) ~ L^(-ln|V'lVlnL)-<i, (1) 

In terms of the variable a = — ln|V'P/lnL, the PDF is 
7^l(q;) ~ where f{a) is the multifractal spectrum, 

i.e. the fractal dimensions of the different a-sets made up of 
the points where j-^ip = L^". The standard f{a) spectrum 



is usually constructed by a Legendre transformation ifiell of 
the scaling exponents T{q) for the generalized inverse partici- 
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pation ratios L'^dV'jp?) ~ L'^W |2l |3l |4l |9l IJI 

The relation (HJ between Vhia) and f{a) suggests a com- 
plete characterization of multifractality directly from the PDF. 
The proportionality in ([T]i contains an L-dependence which 
can be naively included as 



(2) 



where ao is the position of the maximum of the multifractal 
spectrum, /(ao) = d. Furthermore, since f{a) < d for all 
a ^ ao, we see that ao is in fact also the position of the 
maximum of the PDF itself. Hence VL{cto) corresponds to 
the maximum value of the distribution, and ao can be eas- 
ily obtained numerically. As for /(a), the value of ao ob- 
tained from the PDF must be L-invariant at criticality, as we 
show in Fig.[T] The estimation for ao = 4.024 ± 0.022 from 
the PDF is in agreement with that obtained from gIPR scaling 
lUl], ao e [4.024,4.030]. From the normalization condition 
we find 7'L(a^]^= L^(")-''da) "\ Using the saddle 



point method lllSll . justified in the limit of large L, we com- 
pute 7^L(ao) ~ Vln L, which holds very well even for small 
L as shown in Fig. [TJa). From the PDF for fixed L the multi- 
fractal spectrum is hence straightforwardly obtained from (|2]i 
as /(a) = d + ln['PL(a)/'PL(ao)]/lnL. Alternatively, if 
/(a) is known, the PDF can be easily generated. We find ex- 
cellent agreement between the singularity spectrum obtained 
from the PDF for L = 30, . . . , 240 and the one obtained from 
the more involved box-size scaling of the gIPR It must be 
emphasized that Vl{c() is always system-size dependent, and 
it is through the /(a) spectrum that all PDFs for different L 
collapse onto the same function, cp. Fig.fllb). Thus the /(a) 
can also be understood as the natural scale-invariant distribu- 
tion at criticality. 

In order to minimize finite-size effects, we can also deter- 
mine /(a) from the PDF using system-size scaling. We note 
that for a given L the number of points per wavefunction with 
a€ [a-Aa/2,a + Aa/2] is7VL(a) = L'^Vl {a) Aa. Hence 
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a=-lnm /InL 



FIG. 1: (color online) PDF at criticality for Aa — 0.04. The gray 
lines correspond to L from 30 (bottom) to 90 (top). Standard devia- 
tions are within the line width. For L 100 and L > 100 we aver- 
age over 2.5 x 10* and 100 states respectively. The vertical dashed 
line marks the mean value for uq — 4.024 ± 0.022 using L from 50 
to 200. The shaded vertical region corresponds to its 95% c.i. Inset 
(a) shows the maximum values Pl(qo) vs L. Standard deviations 
are contained within symbol size. The solid line is the fit a \n{L/l)'', 
with a = 0.298 ± 0.002, b = 0.489 ± 0.006. Inset (b) shows the 
collapse of all the PDF from L = 30 to 240 onto the f {a). 



the following normalized volume of the a-set J\fL{a) 
i'^PL(a)/7'L(ao) obeying 
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can be used to exttact /(a) from a series of systems with dif- 
ferent L. In Fig. [2] we compare the multifractal spectrum ob- 
tained using PDF scaling ^ with the one from gIPR scaling 
lUSl for L e [20, 100] and having 2.5 x lO'' critical states for 
each size. We find very good agreement between both, as well 
as between the numerical PDFs and those generated from the 
gIPR/(a) [Fig.ma)]. 

Numerically, the PDF is approximated by the histogram 

VL{a) = (6l(Aa/2- |a + ln|V'»lVlni|))/Aa, (4) 

Aa — >0 

where 9 is the Heaviside step function and (.) involves an 
average over the volume of the system and all realizations 
of disorder. To minimize the uncertainty in the small 
which becomes greatly enhanced in terms of ad^/jip), the 
amplitudes used for the histogram are those obtained from a 
coarse-graining procedure of the state using boxes of linear 
size I — 5. Therefore the system size in all equations is the 
effective system size L/l. The uncertainty of the PDF value 
is estimated from the usual standard deviation for a counting 
process as cr-Pi(a) = y^VL{ct)/NwL'^Aa, where is the 
total number of states in the average, and for the a values a 
constant ctq = Aa/3 is assigned. We note that this proce- 
dure assumes uncorrected a's and hence ji/'P's; this is only 
true between different disorder realizations, but not necessar- 
ily within each state. Hence the errors of the PDF are proba- 
bly somewhat underestimated and the small uncertainty of the 
/(a)-values obtained from PDF scaling in Fig.[2]must be in- 
terpreted carefully. There also exists another source of error 
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FIG. 2: (color online) Comparison between f{a) obtained from 
system-size scaling of the PDF for L £ [20, 100] with 2.5 x 10* 
states for each L and the spectrum obtained from the same set of 
data using system-size scaling of the gIPR. Only one every second 
symbol is shown for clarity. Standard deviations are within symbol 
size when not shown. The inset shows the numerical PDF (black) 
and the PDF calculated from the multifractal spectrum obtained from 
gIPR scaling (gray). The bottom panel gives the values of the linear 
correlation coefficient and quality-of-fit parameter Q for the /(q) 
obtained from the log-log linear fits of Eq. ([3]l. 



difficult to quantify, namely, how much the histogram for fi- 
nite Aa deviates from the real PDF when Aa 0. In spite 
of this, the PDF method is easy to implement numerically and 
hence a valid alternative to the more demanding gIPR scaling 
techniques. 

Symmetry relation for the PDF. The symmetry relation, 
f{2d — a) = /(a) + d — a iQ] for L oo, implies the 
existence of a symmetry also for the PDF which should hold 
for large enough system sizes. 



a] 



(5) 



In terms of the wavefunction amphtudes, Eq. ^ reads 

Vl{L-^'^/\iIj\^) = {L'^\ij\'^frLi\M^)- The latter relation 
establishes that at criticality the distribution in the interval 



in the region L^' 



L is indeed determined by the PDF 
^ \ip\'^ < 1. We carry out a nu- 



merical check of the symmetry relation ^ by evaluating 
SVLia) = VLia) - L°'-'^rL(2d- a) accounting for the dis- 
tance between the original PDF and its symmetry-transformed 
counterpart at every a, as well as the cumulative difference 
S{L) = J^"^ da\6VLict)\- The symmetry-transformed PDF 
for L = 100 and the evolution of 6Vl{c() for different L are 
shown in Fig. [3] We find that the symmetry relation is better 
satisfied as L increases [l^j, and the improvement can be 
roughly quantified as S{L) ~ /^-0-537 shown in Fig.[3jb). 

Non-gaussian nature of the PDF. The parabolic ap- 
proximation for /(a) in d = 2 + e, 11911 implies a 
gaussian approximation (GA) for the PDF, Vj^^{a) — 
^yi^^j^y4^^-[a-(d+e)]V4£^ At first glance, the PDFs in 
Fig. [U might indeed appear roughly gaussian and in Fig. |4] 
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FIG. 3: Degree of symmetry of the PDF using Aa — 0.04 for differ- 
ent system sizes with 2.5 x 10^ critical wavefunctions for each size. 
The shaded region in the background marks the 66% c.i. of the val- 
ues for L = 100. The lines show the evolution in average for each L. 
The inset (a) displays the symmetry transformed PDF for L = 100 
and its 95% c.i. (shaded region in the background). The inset (b) 
shows a log-log plot of the cumulative difference S{L) normalized 
to its maximum value (minimum L). The solid line corresponds to a 
linear fit with slope b = -0.537 ± 0.016. 

we show gaussian fits of the PDF for L = 100, obtained via 
a usual minimization taking into account the uncertainties 
of the PDF values. However, the quality-of-fit parameter Q, 
which gives an indication on the reliability of the fit, is ridicu- 
lously small {Q < 10-160000-) since the individual standard 
deviations of the PDF values may have been slightly underes- 
timated we also study how Q behaves when we intentionally 
increase the error bars by a factor n. As shown in Fig. Uc) 
one would have to go to unreasonable high values of ti ^ 60 
to accept a gaussian nature for the PDF as plausible. The de- 
viation from Vj^^{a) is also noticeable. Hence our statistical 
analysis confirms that the PDF is non-gaussian in agreement 
with the observed non-parabolic nature of f{a) at the 3D MIT 
SIB. 

Rare events and their negative fractal dimensions. The 
volume of the a-set, AfL{ct) for a given L, gives the num- 
ber of points in the wavefunction with amplitudes in the range 
\^pi\'^ e [2.-"-^"/2,L-"+^"/2]. It scales with the system 
size as Mhioi) ^ Vln L L^''°'K The negative values of f{a) 
ifioll correspond then to those a-sets whose volume decreases 
with L for large enough L. Physically, the negative fractal di- 
mensions at small a are caused by the so-called rare events 
containing localized-like regions of anomalously high 
at criticality. The probability of finding them likewise de- 
creases with L. In Fig. [5] we show examples of rare eigen- 
states. Due to the finite size term -y/ln L, the threshold 
[where f{a-) = 0], below which the decreasing behaviour 
of Ml (a) with L is detected, will change with the system size 
itself, and so the normalized volume (|3]l of the a-set must be 
used. In Fig. |6]we show the behaviour of Ml (a) vs L for two 



FIG. 4: (color online) Gaussian fit (y^a/^ e'"'^'''''^') of the PDF 
obtained from 2.5 x 10'* critical states of size = 100^ us- 
ing Aa — 0.04. The dashed-dot line corresponds to the fit with- 
out taking into account the uncertainty of the PDF values (a = 
0.815 ± 0.004, b = 3.9731 ± 0.0021). The dashed line shows the 
fit obtained when the uncertainty of the PDF points is considered 
(a = 0.840 ± 0.005, b = 3.942 ± 0.004, Q < iQ-i""""). The 
solid gray line shows the GA for e — 1. Insets (a) and (b) are blow- 
ups of the maximum and the right tail of the PDF. Inset (c) shows the 
quality-of-fit parameter Q when the standard deviation of the PDF- 
points is multiplied by the factor n. 




FIG. 5: (color online) Rare eigenstates for the 3D Anderson model 
at S = and Wc = 16.5 for L = 100. The sites with probability 
I t/jj P larger than L ^ are shown as boxes with volume | ?/)| | L'^ . The 
grayscale distinguishes between different slices of the system along 
the axis into the page. The biggest boxes with black edges enclose 
the site with the maximum normalized amplitude: Itpil^ — 0.4484 
(left) and 0.3617 (right). 

values of a corresponding to a positive and a negative fractal 
dimension. The exponential decreasing of the volume of the 
a-set for /(a) < is clearly observed from the PDF values. 
In Fig.|6ja) it is demonstrated how the normalized volume of 
the a-set becomes scale invariant at a_ and thus Ml (a) = 1. 
The opposite tendencies with L at each side of a_ can also be 
seen. The estimated value for a_ e [0.644, 0.677] from the 
PDF scaling agrees with the result obtained using the multi- 
fractal spectrum from gIPR scaling, a_ = 0.626±0.028lfl0ll. 

Termination points in /(a). The fate of the /(a) spec- 
trum at a = and a = 2c? is currently under debate in 
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FIG. 6: (color online) NL{a) vs L for a = 0.45 (filled squares, 
left y axis) and 0.75 (open squares, right y axis). The solid lines 
are the curves where f{a) is obtained from the log-log lin- 

ear fits of Eq. The shaded regions marlc the 66% c.i. of the 
^/('^) values. The inset shows ln[NL{ce)] vs a for different L in 
the range [20 (black circle), 100 (black triangle)]. Standard de- 
viations are within symbol size, whenever not shown. The verti- 
cal dashed line corresponds to the mean value of the crossing point 
a- = 0.661 ± 0.008 and the shaded bar marks its 95% c.i.. 

the literature lfloil20ll due to the emergence of singularities at 
these points. At present, it is not clear whether /(a) contin- 
ues towards — oo or terminates with finite values. The phys- 
ical consequences of the absence of termination points (TP) 
are: (i) jV'iP > L^^"^ at criticality since 2d would be an up- 
per bound for a, and (ii) the probability to find the most rare 
event, namely the most extremely localized state (iV'ioP ~ 1' 
corresponding to a = 0), at the critical point must always be 
zero independently of the system size [7^1,(0) = 0]. In princi- 
ple the PDF can be used to look for TPs both at a = and 2d. 
However, a reliable analysis in the vicinity of a = requires 
a huge number of disorder realizations; relying on the sym- 
metry relation 101 a study around a ^ 2d is more appropriate. 
For 1 a < 2d and as long as there is a TP, the PDF admits 
the series expansion 

VLia) ~ VLiao)Lf^^'''^-'' [1 + q2d{a - 2d)lnL] , (6) 

where q2d = f'{c()\a=2d- The existence of a TP requires 
,f{2d) and q2d to be L-independent, for large enough L. In 
Fig.|7]we show the values of f{2d) and q2d obtained for dif- 
ferent L. The value of f{2d) seems to reach a saturation for 
large L although the numerical analysis cannot exclude a very 
slow decreasing tendency. A similar result is found for q2d- 
Still larger L and more states are needed to decide the fate of 
f{a) at and 2d. 

In conclusion, we have shown here that a PDF-based study 
can be a valid companion approach to the standard multifrac- 
tal analysis, giving complementary and new information about 
the critical properties of waves at Anderson-type transitions as 
well as offering a conceptually simpler viewpoint. 
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FIG. 7: Values of /(2d) (filled squares, left y axis) and q2d = 
f' {a)\a=2d (open squares, right y axis) as functions of system size. 
The inset shows the fits of Eq. ^ for different L from 20 (top) to 
100 (bottom). Standard deviations are within symbol size. 
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